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POISSON VORONOI APPROXIMATION 

By Matthias Heveling^ and Matthias Reitzner 

Vienna University of Technology 

Let X be a Poisson point process and 7^ C K'' a measurable set. 
Construct the Voronoi cells of all points x £ X with respect to X, and 
denote by vx{K) the union of all Voronoi cells with nucleus in K. 
For K a compact convex set the expectation of the volume difference 
V{vx{K)) - V{K) and the symmetric difference V{vx{K)l\K) is 
computed. Precise estimates for the variance of both quantities are 
obtained which follow from a new jackknife inequality for the variance 
of functionals of a Poisson point process. Concentration inequalities 
for both quantities are proved using Azuma's inequality. 

1. Introduction. Let X be a stationary Poisson point process of intensity 
A. Denote by vx{x) tlie (random) Voronoi cell of x with respect to the point 
set X U {x}, that is, 

vx{x) = {z ^W^: \\z — x\\ < \\z — y\\ for all y G X}. 

We call X the nucleus of the Voronoi cell vx{x). The set of all Voronoi cells 
vx{x)-, X ^ X, is the Poisson- Voronoi tessellation of W^. For a given set 
A C M"' we consider the Poisson- Voronoi approximation vx{A) of A which 
consists of all Voronoi cells with nucleus in A, 

vxiA) := IJ vx{x). 

The set vxiA) is a random approximation of A. In this paper, we discuss 
the quality of this approximation for a convex set A. In particular, we are 
interested in the convergence of vx (A) to A when the intensity of the Poisson 
process tends to infinity. 
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The general problem whether vxiA) approximates A for "complicated" 
sets A seems to be difficult; only partial answers are available (see Khmal- 
adze and Tor onj adze [7] and Penrose [14]). 

Here we concentrate on the case of a compact convex set K with nonempty 
interior, and its approximation vxiK), where we derive precise estimates for 
the volume V{vxiK)) and the volume of the symmetric difference of K and 
vx{K). 

Theorem 1. Let X be a stationary Poisson point process of intensity 
X. If K is a convex set of volume V{K) and surface area S{K), then 

EV{vx{K)) = V{K) 

and 

YaTV{vx{K)) < Cd\~^~^''^S{K) 

with a constant Cd independent of K and A. Further there are constants 
c{K),A{K) such that 

F{\V{vx{K)) - V{K)\ > t^\~^-^/dS{K)) 

< ciK)e-<''^''-^^''''"' + 16VdX-''+'SiK) 
with c'^ = 2-43-2'='d-'^-^/2 A > A{K) and k>2. 

Our proof of the theorem relies on the stationarity of the process. Its first 
part can be generalized to nonstationary Poisson processes with an absolute 
continuous intensity measure with density function Xf{x) > (with respect 
to Lebesgue measure). In that case the volume on the right-hand side of 
the formula is replaced by the integral of the density function / of the 
intensity measure over K. Observe that if this density vanishes on a set of 
positive measure, then we do not even have KV{vx{K)) — > y{K) for certain 
sets K. The second part of the theorem concerning the variance and the 
tail probability is also subject to generalization and holds for nonstationary 
Poisson processes with bounded density function. The present form of the 
theorem, however, gives fastest access to our asymptotic results. 

Since the expectation of the volume of the Poisson- Voronoi approximation 
coincides with the volume of the convex set itself, it is natural to ask for the 
volume of the symmetric difference KAvx (K) = {K \ vx {K)) U {vx (K) \ K) ■ 
It is known that for any bounded Borel set yl C M'^, one has 

V{AAvx{A)) ^0 

almost surely as A ^ oo. This was proved for d = 1 in Khmaladze and Toron- 
jadze [7], and by Einmahl and Khmaladze [3] for any bounded Borel set 



POISSON-VORONOI APPROXIMATION 



3 



AcR'^ with V{Ae) V{A) for e ^ 0, where Ae = A + eB'^ is the Minkowski 
sum of A and the bah eB'^. The general case was proved by Penrose [14]. 

Here we concentrate on rates of convergence for convex sets A. We derive 
an asymptotically precise formula for the expectation and estimates for the 
variance and the tails. We believe that the estimates for the variance in 
Theorems 1 and 2 are best possible up to the choice of the constant. Denote 
by Kd the volume of the unit ball in W^, and by r(-) the Gamma function. 

Theorem 2. Let X he a stationary Poisson point process of intensity 
A. If K is a convex set of volume V{K) and surface area S{K), then 

(1) ¥y{K^vx{K)) = ceA-i/'^5(k)(i + 0{\~^'^)) 

with CE = |K^^^'^Kd_ir(^). And 

YaTV{KAvx{K)) < CdX'^-^/'^S{K) 

with a constant c^ independent of K and A. Further there are constants 
c{K),A{K) such that 

F{\V{KAvx{K)) - EV{KAvx{K))\ > t^ \-^-^/'^S{K)) 

< c{K)e~'^'^''^^'^^^'"'"" + 16VdA-^+i5(K) 

with = 2-43^2d^-<i-l/2 ;y > ^ > 2. 

Remark 1. A precise estimate for the error term in (1) is given in 
Section 5, Theorem 6. 

Remark 2. In both theorems the estimates for the tail probabilities are 
stated for \ >K{K) sufficiently large. Set r = Ay/d{k\-^\n\Y/'^ . Then A is 
sufficiently large, if r < 4 and the volume of the parallel set V{dK + rB'^) is 
bounded by ArS{K). 

Remark 3. Since in both cases, the results for the expectation and 
variance only depend on the volume and surface area of K, they possibly 
hold for more general classes of sets. Yet our methods of proof make essential 
use of the convexity of K (in particular Hadwiger's characterization theorem; 
see Section 2). In view of applications it would be of interest to extend our 
results to sets in the convex ring. 

Remark 4. Jeulin posed the following interesting problem: Is it better 
to approximate V{K) by a single realization vx{K) where X is a Poisson 
point process of intensity A = /cAq, or by the mean value of k realizations of 
vXiiK) where Xi, i = 1, . . . ,k, are independent Poisson point processes of 
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intensity Aq? Both estimates are unbiased. Yet the approximation by a single 



Our results have applications in nonparametric statistics (see Einmahl and 
Khmaladze [3], Section 3) and image analysis (reconstructing an image from 
its intersection with a Poisson point process; see [7]). Also the connection to 
quantization problems is obvious; see, for example, Chapter 9 in the book of 
Graf and Luschgy [4] which gives an excellent introduction into this topic. 
Quantization problems are connected to problems of numerical integration: 
approximate the volume Vd{A) = dx of a set A using its intersection 
with a point process X r\ A. As shown in Theorem 1 the volume of the 
Poisson- Voronoi approximation vx{A) is an unbiased estimator for Vrf(A), 
even for Borel sets A which is pointed out in Section 4. An estimate for the 
volume of A is also obtained if the number of points X{A) that fall into A 
is counted instead of the volume of the Poisson- Voronoi approximation. By 
the definition of X one has 



When the variation coefficient is regarded as a measure for the quality of the 
respective volume estimators, then using Poisson- Voronoi approximation is 
more precise than counting points of the Poisson point process — at least for 
convex sets A and stationary Poisson point processes. It would be of interest 
to show that this is a general principle even for arbitrary point processes, 
for example for random lattices X. 

An interesting open problem is to measure the quality of approximation 
of K by vx{K) using the Hausdorff distance between both sets. We are not 
aware of any results in this direction. 

For basic facts from integral geometry, stochastic geometry and Voronoi 
tessellations which are not explained in the following, we refer the reader to 
[13, 16, 17, 18]. The employed notions and results from the theory of convex 
bodies are found in [15]. 

We work in d-dimensional Euclidean vector space W^, with norm d{x, y) = 
\\x — y\\, and for closed sets K <ZW^ distance d{x, K) = min(d(2;, y)\ y ^ K). 
Its unit ball, {x G M'^:|2;|| < 1}, is denoted by B'^, and S'^^^ is the unit 
sphere. The space of convex bodies (nonempty, compact, convex subsets) 
in is denoted by IC^ and the space of locally finite point sets in is 
denoted by N. 

For a stationary Poisson point process, as usual, X denotes the simple 
counting measure as well as its support, that is, X{A) and card(Xn^) have 
the same meaning. Its intensity measure = EX(-) (E denotes mathemati- 
cal expectation) is given by 



estimate should be preferred, since \ai{vx{I^ 
whereas Var(^ X]i=i "^x, (-f^)) is of order k~^\ 




¥.X{A) =\ai X{A) = XV{A). 
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2. Valuations and Delaunay triangulations. A major tool for proving our 
results is valuations. A functional /i : Kf^ ^ M is called a valuation, if for every 
K,LelC'^ with KULgIC^, 

li{K U L) + n L) = ^l{K) + 

holds. A monotone valuation satisfies ^{K) < fJ-{L) if K C L. Valuations play 
an important role in convex geometry and integral geometry; for further 
references see [9, 11] and [12]. One of the most important results in this field 
is the following characterization theorem by Hadwiger [5]: 

A functional fi : IC^ ^ M is a monotone and rigid motion-invariant valua- 
tion if and only if there are constants cq, ci, . . . ,Crf {uniquely determined by 
fi) such that 

fi{K) = CdVd{K) + --- + coVo{K) 

for every K G /C^ . 

Here Vq{K) , . . . ,Vd{K) are the quermassintegrals of K. In particular, 
Vd{K) equals the volume V{K), 2Vd-i{K) is the surface area S{K), and 
Vq{K) is a multiple of the Euler characteristic. For a modern proof of this 
theorem, see Klain [8]. 

In the following sections a valuation depending on the Delaunay mosaic 
of X turns out to be of importance. To this end denote by £x the edges of 
the Delaunay mosaic of X , that is, those segments [x,y] with x,y £ X such 
that vxix) Civxiy) is a facet of vx{x) and vx{y)- Set 

nK[x,y] = 21{[x,y]nK^ 0) - l{x e K) - l{y G K), 

where l(-) denotes the indicator function. Observe that for points x,y £ 
X, with probability 1, nK[x,y] is the number of connected components of 
[x,y] \K if [x,y] meets K, and nK[x,y] = otherwise. 

Theorem 3. Let / : N x M'' ^ [0, oo) be measurable. Define a functional 
fi:IC^^M 62//i(K):=EE[.,,]6^^-(/(^,x) + /(X,y))nA'[a;,y]. Then 

(2) fi{K) = CfiX)SiK), 

where S{K) is the surface area of K. Moreover, if f has the scaling property 
f{tip,tx) = t°'f(ip,x), then there is a constant Cf such that 

Proof. First we will show that fj, can be expressed as the difference of 
two auxiliary monotone valuations. We define 

ui{K):=E J2 mx,y]nK^0){f{X,x)+f{X,y)). 
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Elementary considerations yield that for any two convex bodies K and L 
such that Kyj L \s convex we have 

l{[x,y\r\K^0) + l{[x,y\r\L^0) 

= y] n (K U L) / 0) + y] n {K nL)^ 0). 

Hence z^i is a valuation and it follows directly from the definition that i^i is 
also monotone and rigid motion- invariant. We define a second functional 1^2 
by 

MK):=^ E m^^K) + l{yeK)){fiX,x) + f{X,y)). 

As above for z^i it is easily seen that z/2 is also a monotone and rigid motion- 
invariant valuation. Then Hadwiger's theorem yields that both ui and 1^2 can 
be written as a linear combination of the Minkowski functionals. Moreover, 
since /i = 2ui — 1^2, we can write ^{K) = J2i=oCi{\)Vi{K). If the dimension 
of K is less than d — 1, then nK[x,y] = for all [x,y] G £x with probability 
1 and thus cq, . . . , Cd-2 are vanishing. Hence 

fi{K) = Cd-i{X)Vd-i{K)+Cd{X)Vd{K). 

If K is of dimension d — I, then IJ-{K) = 2ui{K) = Cd-i{\)Vd-i{K) . Sup- 
pose that P is a polytope with facets F G J'{P) and nonempty interior. Then 
the valuation //(•) can be written as 

^Ji{P)= E E {f{X,x) + f{X,y)) 

= Cd-lWVd-l{P) 

and thus = which proves (2). 

For the second claim of the theorem the scaling property of the Pois- 
son process is used. Write fi\{K) = EJ2ix,y]eex(f(-^'^) + f{X,y))nK[x,y] 
to emphasize the dependence on the intensity A of the point process X . For 
t>0 replace x,y by x/t,y/t. Then 

^lxiK)=E E t-''ifitX,x) + f{tX,y))ntK[x,y] = t-''fi,-,^{tK), 

[x,y]&£tx 

and, together with (2), we obtain 
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3. A jackknife estimate for the variance of functionals of a Poisson point 
process. To get an estimate for the variance of a function S{X) we rewrite 
the Efron-Stein jackknife inequality [2] (see also Efron [1] and Hall [6]). In 
the form we are interested in this is possible if there are no far-reaching 
dependencies. This is made precise in the following assumptions: 

For a locally finite subset 1" S M"^ we call R{Y) the radius of influence 
of a function S{Y), if there is a function /:R — > M such that for arbitrary 
locally finite point sets DcR'^ and D° C (M"' \ R{Y)B'^), we have: 

(Al) S{Y)=S{YUD°); 

(A2) \SiY)-SiYuD)\<f{RiY)). 

Hence the influence on S{-) of additional points can be estimated by the 
function f{R{Y)), and additional points outside R{Y)B'^ are negligible. This 
notion is close to the notion of stabilization used in previous work, for ex- 
ample, in Penrose [14], where stabilization at the origin refers to condition 
(Al) whereas here we need in addition bounds on the costs of adding points 
close to the origin. 

We are interested in the case when Y is the realization of a Poisson point 
process X and the moments of f{R{X)) are bounded. 

Theorem 4. Let X (resp. X~^), be a Poisson point process of intensity 
X, resp. A(l -|- ^). Let S':N — > M be a measurable function on the space of 
locally finite point sets in M*^, and let R{X) be a radius of influence of S{X). 
If¥.{f{R{X)fR{Xf'^) exists, then 

Var5(X) < Jirn^E ^ {S{X+ \{x}) - S{X+)f . 

Proof. We start with recalling the Efron-Stein jackknife inequality in 
its usual form. Let Yi be independent identically distributed random vari- 
ables defined on some probability space, i = 1, . . . , m -|- 1. We write y^*^ for 
(Yi, . . .,Yi^i, Fj+i, . . .,Ym+i). If S'(yi, . . .,Ym) is any real symmetric func- 
tion of m random variables, an estimate for the expectation of S{-) is given 
by 

1 m+l 

s=——y s{Y^^). 

m + l ^ ^ ' 

1=1 

The Efron-Stein jackknife inequality then says that the natural estimate for 
the variance '}2{S{Y^^'^) — S)"^ overestimates the real variance: 

m+l 

Var5<E^(5(yW)-5)l 
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Since the right-hand side increases if we replace the mean S by any other 
function T = T(Yi, . . . , Ym+i), we also have 

m+l 

(3) YaiS<Ej2iS{Y^''^)-Tf- 

i=l 

Let Xi, . . . , Xm+i be independent Poisson point processes in M*^ of in- 
tensity A/m, set X« = Uj=i,...,i-M+i,...,m+i X = X('^+'^\ and X+ = 
Uj=L^ Xj which are Poisson point processes in M*^ of intensity A, or A(l + ^) 
respectively. Since by assumption S{Xi, . . . ,Xm) = S{Xi U • • • U X^) is a 
symmetric function in the Xi, the Efron-Stein jackknife inequality (3) with 
T = S{X+) tells us that 

m+l 

Var S{X) <EJ2 - SiX+)f 

i=l 
m+l 

= E ^ eW(S(xW) - 

i=l 

where E(*)(-) abbreviates E(-|X(*)). 

For the next step fix i and denote the radius of influence of S{X^^^) by 
R{X^^^) = R. Observe that this implies that R is independent of X^. So we 
may apply conditions (Al) and (A2) with Y = X'^'\ D° = Xi\RB'^, D = Xi. 

If m is large, then with high probability at most one point x £ Xi is in 
RB'^ and thus may have influence on S{X~^). This is made precise in the 
following. We decompose the expectation according to the value of Xi{RB'^): 

eW(S(X»)-S(X+))2 

oo 

= Y,^^'H{S{X^''^) - S{X+)fl{Xi{RB'^) =n)). 

n=Q 

For Xi{RB'^) G {0, 1} we use (Al) with D° = Xi\ RB'^ and obtain 
X:E»( E {S{X+\{x])-S{X+)fl{X,{RB^) = n)\ 

< g eW I Y: {S{X+ \ {x}) - S{X+) fliX,{RB') = n)] 

n=0 \x€Xi J 

= E« E ('^(^^ \ {^}) - S{X+)f. 

x£Xi 
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For Xi{RB'^) > 2 we use (A2) with D = Xi, the estimate X;n=2 ~re~^'<^, 
and obtain 



n=2 n=2 

\2m^ UT}d\2 



£ E»((S(XW) - S{X+)fl{X,{RB'') = n)) < £ f {Rf^MRB'') = n) 

2X^ V(RB<^f 

since the intensity of Xi equals ^ . 

Combining our results, summing over z = l,...,m + l, and using that the 
radii of influence R^X^"^^) are identically distributed, gives 



(4)Var5(X)<E V {S{X+ \{x]) - S{X+)f + —¥.{f{RfV{RB^^'^' 
^-^ , m 
X&X+ 

and thus proves the theorem. □ 



In the next section we use Theorem 4 for functionals with moments con- 
tinuous in the intensity of X. From (4) we obtain in this case the following 
corollary. 



Corollary 5. Let X he a Poisson point process. Let SiN^R he a 
measurable function on the space of locally finite point sets in M*^, and let 
R{X) he a radius of influence of S{X). IfE,S(X), ES'^(X) are continuous 
in X, and if W.{f {R{X)f R{Xf'^) exists, then 

Var5(X) < E ^ {S{X \ {x}) - S{X)f. 



We want to remark that the Slivnyak-Mecke formula for a Poisson point 
process allows to rewrite our theorem in the following way: 

Var S(X) < A / lim E(S(X+) - S(X+ U {x})f dx. 

We conjecture that the following more general theorem holds: 



Conjecture. Let X be a Poisson point process. For any measurable 
function 5" : N — > M on the space of locally finite point sets in we have 



\avS{X) < E ^ {S{X \ {x}) - S{X)f. 
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4. Volume difference. In this section we are interested in the difference of 
the volume of vx{K) and K. We state the mean value and prove an estimate 
for the variance. The large deviation inequality is proved in Section 6. 

It can easily be shown that 

(5) mr{vx{K)) = v{K). 

This follows either from Campbell's theorem (see, e.g., the book by Schneider 
and Weil [17], page 128) or using Hadwiger's characterization theorem and 
an argument similar to that of Theorem 3. Formula (5) holds for all Borel 
sets without any convexity assumptions. 

To get a bound on the variance of V{vx{K)) we use the Efron-Stein 
jackknife inequality in Corollary 5. This states that 

\s.rV{vx{K)) < E ^ {V{vx\{^,}{K)) - V{vx{K))f 

x&X 

if for some radius R{X) of influence the moment E(/(i?(X))^i?(X)^"') ex- 
ists. (Observe that the moments of the functional we are interested in are 
continuous in A.) 

Thus we have to estimate the volume of those Voronoi cells with centers 
X S X, which partly may change from exterior points to interior points or 
vice versa if x is removed. Assume that X n dK is empty which happens 
with probability 1. 

If for X G X n A' all neighbors of the Voronoi cell vx{x) are also contained 
in K, that is, if for all [x,y\ G £x we have y £ K, then Vx\{x}(,K) = vx{K). 
The same argument applies if the point x and all its neighbors are outside 
K. Hence of interest are those points x £ X such that there exists an edge 
[x, y\ S £x with [x, y] H dK ^ in which case 

\V{vx\{x}{K)) - V{vx{K))\ < V{vx{x)). 

Defining nj<-[x,y] as in Section 2 and noting that nft-[x,y] > l{[x,y] H dK ^ 
0) we thus see that 

YaTV{vx{K))<E J2 nK[x,y]{V{vx{x)f + V{vx{y)f). 

[x,y]&£x 

By Theorem 3 with a = 2d we immediately obtain 

\aTV{vx{K)) < CdX'^'^^'^SiK) 

which is the variance estimate of Theorem 1 . 

It remains to define the radius of influence R{X), and to show that 
K{f{R{X))'^R{Xf'^) exists. Define a (random) number R' = R'{X) as the 
smallest number fulfilling 

(6) \Jvx{z)cR'B^ 
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and let R{X) = 3R' . [Recall that vx{K) is the union of all Voronoi cells 
vxix) with nucleus x £ X CiK, whereas here we estimate the influence of all 
zeK.] 

As for assumption (Al) we have to show that any point set D° which does 
not meet 3R'B'^ has no influence. Indeed, if vxiK) ^ vxuD°{K), then there 
are points x G X Ci K and y £ D° with [x,y] G £xuD°- Hence the Voronoi 
cells of x and y would have points in common. This is impossible since by 
the definition of R' we have vx{x) C R'B'^, but the midplane between x and 
y does not meet the ball R'B'^. 

As for assumption (A2), it follows from (6) that for any point set D 

vxud{K) C vxu{DnK){K) c R'B''. 

Thus the difference between vx{K) and vxuDiK) is bounded by the volume 
of R'B'^ and assumption (A2) is fulfilled with fiR{X)) = V{R'B'^). 

Finally we have to show that E(ii'^'^) is finite. Denote by Rk the smallest 
radius such that K C RkB'^. By definition, if i?' > r for some r > Rk + 
then there is a point y € rS"^'^ with d{y, dK) < d{y, X), that is, X{B{y, [r — 
Rk))) = 0. We cover the ball rB'^ by disjoint cubes Ci of sidelength 1 
with center Zi and obtain 

nR' > < E ^i^y e n ■ ^(^(2/' - ^k))) = 0) 
1=1 

2dj.d 

< E G Ci : X{B{zi, {r-Vd- Rk))) = 0) 
1=1 

Thus all moments of R' exist and are finite. 

5. Symmetric difference metric. In this section we investigate the vol- 
ume of the symmetric difference of vx{K) and K, 

V{KAvx{K)) = V{K \ vx{K)) + V{vx{K) \ K). 

We determine the expectation (Theorem 6), and prove an estimate for the 
variance. The large deviation inequality is proved in Section 6. 
First we show that 

V{K/\vx{K)) = c'a_,SiK)X-^/'' + o(A-i/'^). 

We start with the volume of vxiK) \ K. The Slivnyak-Mecke formula gives 
for X G M"^ \ K 



P(x e vx{K)) = ¥{3y eXnK:xe vxiy)) 
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y&XnK 

(7) 

= X[ e-AV(B(x,d(x,s/)))^y 



since x G Vxu{y}iy) if the intersection of X with the bah of radius d{x,y) 
centered at x is empty. Precisely the same argument shows that for x £ K 

(8) F{x^vx{K))=E y l{xevx(.y)) = X f e-^^^^^^'-^^^'^))) dy. 

yeX\K 

Combining (7) and (8) we obtain 

KV{KAvx{K)) =E [ l{x£ KAvx{K)) dx 

= [ F{x£vx{K))dx+ [ F{x^vx{K))dx 

= 2X [ [ e^^^^ll^^^ll^dydx. 
Jk''\k Jk 

We use the Blaschke-Petkantschin formula (see, e.g., [16]) which transforms 
the integration of the tuple {x,y) with respect to Lebesgue measure into 
integration of {x,y) with respect to the (one-dimensional) Lebesgue measure 
on the line E which is the affine hull of the two points, and then integrate 
with respect to the set £f of all lines in using the normalized Haar 
measure v on the set of all lines: 

EV{KAvx{K)) = dXKd [ I I e-^''''\\y-''\\''\\y - x\\'^-'^ dydxdv{E). 

Jef Je\k Jedk 

Identify E with M and EnK with the interval [0, 1] of length l = l{E).U 
I > 0, we obtain for the inner integrations 

^~^^d\^-y\''\x-y\d-Uydx = 2 / c'^'^''^^-^^" (x - y)^"^ dx dy 
,1] Jo Jo Jl 

= ^(AACd)-^ fe-^'^'^y'dy 

2 







(AKrf)-i-i/'^rQ)(i-5(A,ii;)), 



2 

d2 
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where 

since ^ — 1 < 0. Thus we have 

EViKAvxiK)) = ^(AKrfj-i/'^rQ) 1{E n A' / 0)(1 - <5(A,i?)) du{E). 



For the main term we obtain by Cauchy's surface area formula 
2 



|(AKrf)-i/'^r i(i^ n / 0) 



= ^(AK^)^'/'A^rf„irQ)5(A). 

To estimate the error term j 5{\,E) dv{E) assume that the origin of the 
coordinate system is chosen in such a way that 

r{K)B'^ C K, 

where r{K) is the inradius of K. Parametrize the hne E hy E = tu + y, t £M, 
where u G S'^~^ is the direction of E and y . The measure v decomposes 
into the uniform distribution u on the sphere S'^~^ , and for u G S'^~^ , into 
Lebesgue measure in the hyperplane . If the hne E meets K, then the 
point y is contained in the projection K\^i_ of the set K onto u^: 

[ l{EnK^0)6{X,E)du{E)= [ [ 6{X,E)dydLo{u). 
Jet Js<i-^ Jk\^j_ 

We introduce polar coordinates y = rv, where v is integrated with respect 
to Lebesgue measure a on S"^~^ n n-*-. Denote by p{v) = px\ j_ (v) the radial 
function of K\y^± in direction v. Because is a convex set, we have that 
for fixed u and v the chord length l{rv,u) is a concave function in r which 
vanishes at the boundary of K\^_l. Hence 

l{rv, u) > 1(0, u) (^1 - > 2r{K) (^1 - 

for < r < p{v). This yields 



e-AK.;'^ < / r^"^ e-^^''^''-(''^'^'-''/P^-^^\''-^drda{v) 

K\ ^ J5'*-inu^ Jo 



OQ 



2d JS'^-inu-L 



d-1 

< 

- 2d 



(AK,)-i/''rQ)r(A)-V,_i(ALx). 
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Using Cauchy's surface area formula again gives 

We summarize our results: 

Theorem 6. If K e JC^, then 

EV{KAvx{K)) = ^(AKd)-^/'^Krf-irQ)5(i^)(l - X-'/^A), 

where < A < '^K~^^'^T{^)r{K)-\ 

The same arguments which led to the bound on the variance of V{vx{K)) 
will yield a bound on the variance of V{KAvxiK))- We use again the Efron- 
Stein jackknife inequality proved in Theorem 4 showing that 

YaiV{KAvx{K)) <EY,iy{vx\{x}{K)AK) - V{vx{K)AK)f 

where the radius of influence R{X) is defined precisely as in (6). Hence we 
already know that E{f {R{X)f R{Xf'^) exists. 

As in Section 4, of interest are those points x ^ X for which there exists 
an edge [x, y] E £x with [x, ?/] n dK / 0. In this case 

\V{vx\{.,}{K)AK) - V(vx{K)AK)\ < V{vx{x)). 

Thus we obtain 

YavV{KAvx{K))<E ^ nK[x,y]{V{vx{x))^ + V{vx{y)?) 

which is the variance estimate of Theorem 2. 

6. Large deviation inequalities. In this section we prove the large devi- 
ation inequalities of Theorems 1 and 2. The essential tool is Azuma's in- 
equality, in particular the method of uniformly difference-hounded functions 
used by McDiarmid [10]. 

A function / : Jli x • • • x ^ M is called uniformly difference-hounded hy 
h if the following holds: for any (yi, . . . , y^) G r^i x • • • x and for any k 
and any y'^^Vt^ we have 



|/(?/i,---,yfc,---,ym) -/(yi,---,yfc,---,ym)| <b. 
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Let Yi, . . . ,Ym. be independent random variables with Yk G ^Ik- Set ^ = 
/(Yi, . . . , Ym) with / uniformly difference-bounded by b. Then McDiarmid's 
bounded difference inequality says that, for any t, 

P(|C-1EC| >t) < 2e-2*'/(™f''). 

To define the random variables Yi,...,Ym we need some preparations. 
Dissect M*^ into cubes Ci of diameter 6 having pairwise disjoint interior. 
Define 5 such that 

V{Ci) = d~'^/^5'^ = kX^Hn\ 

with k >2. Assume that the cubes are numbered in such a way that for 
i = 1,. . . , mg the cubes Ci have nonempty intersection with dK + 36 B'^, and 
that Ci is disjoint from dK + 35B'^ for i > m^. Since for i = 1,. . . ,ms the 
cubes Ci are contained in dK + A6B'^, we see that 

mg 

J2 V{Ci) = d~'^/^5'^ms < V{dK + A5B'^). 

i=l 

Note that there is a @{K) G (0, 1] such that for e < @{K) we have 

(9) V{dK + eB'^) < AeS{K) . 
Thus for 4:6 < Q{K) we obtain 

(10) ms<16d'^/^S{K)6-'^+\ 

For A large, in each of the cubes Ci at least one point of the Poisson point 
process is contained with high probability. To make this precise denote by 
A the event that for all i = 1, . . . ,ms we have X{Ci) > 1. Since X nCi is 
Poisson distributed we have 

(11) FiA") = 1 - (1 - e-^^(^»))"^^ < mse-^^^^^^ = m^A'^ 

We assume in the following that each cube Ci, i = 1,. . . ,ms, contains at 
least one point. This implies that ii x € X has distance at most 6 to dK, 
then 

(12) vxix)cBix,6). 

To prove this inclusion observe that for any point y £ Ci, i = l, . . . , ms, the 
distance of y to one of the points X CiCi is at most 6, the diameter of the 
cube Ci. Thus if y is contained in some Voronoi cell vxix) with x G X, then 
X is the nearest point of X to y, and we have 

(13) \\y-x\\<6. 

Now let X G X have distance at most 6 to the boundary of K. Assume 

yeB{x,26)\B{x,6). 
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Then the distance of y to the boundary of K is at most 3(5 and thus y 
is contained in some Q, i = 1, . . . ,7715. By (13) the distance from y to the 
nearest point oi X r\Ci is at most 5 and hence ||y — > (5 imphes that 
y ^ vx{x). This proves (12) since the Voronoi cell vx{x) does not meet 
B{x,26) \ B{x,5) and is connected. 

Since all Voronoi cells meeting the boundary of K have circumradius at 
most 6, they are contained in dK + 25B'^. By (13) the centers of neighbors 
cells, having boundary points y in common with these cells, have distance 
at most 35 to dK and thus are contained in Cj, i= l,...,m5. In other 
words, the set of all Voronoi cells meeting the boundary of K only depends 
on X n Cj, i = 1, . . . , m^, and is independent of all points of X outside the 
cubes Cj, i = 1, . . . , m^. 

For i = 1, . . . , define the random points Yi by 1^ = X n Ci. If C = /(^) 
is a function depending only on those Voronoi cells meeting the boundary 
of K, then C, depends only on Yi, C,= f{X) = f{Yi, . . . , Ymg)- In the cases we 
are interested in, we have either 

C = f{Yi,...,Y^,) = V{vx{K)) - V{K) 
= E E V{vx{x)\K)- E y{vx{x)r^K) 

i=l XxeYinK x&Yi\K 

or 

C = f{Y^,...,Ym,) = V{K^vx{K)) 

ms / 

= E E Vivxix)\K)+ E Vivx{x)nK) 

i=l \xeY^nK xeYi\K 

In both cases it follows from (12) that, replacing Yi by some nonempty finite 
subset Y- cCi, we have 

|/(.. .,y„. ..)-/(.. .,i^',...)l<m + '5i?') 

and thus b = ^'^6'^ is the required difference bound. Now McDiarmid's theo- 
rem tells us that 

P(|C - E(C|^)| > t\A) < 2e-2*V(m.fc^). 
It follows from P(-) < P(-|A) +P(A^) and from (11) that 

P(|C - IE(C|^)| >t)< 2e^2tV("^*fe2) + ^^A-^ 
In the last step we replace E((^|74) by E^. We use the elementary inequality 
|EC - E{C\A)\ < mCHA)) - KiC\A)\ + EiCHA')) 
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< E(C|^)P(^'=) + ^JE{C)F{A'=) 

<iE{C\A) + .jE{e))^/nA^), 

where the second Hne follows from Holder's inequality. Since, conditioning on 
A, all Voronoi cells meeting the boundary of K are contained in dK + 26B'^ , 
we have by (9) 

E{C\A)<85SiK). 

And for E^^ the bounds on the expectation and variance yield immediately 
that 

]EC^<C2(K)A-2/'^ 
for A > 1 (which follows from 6 <1). Thus 

2|EC - E(C|A)|2 mj^b'^ < cz{K)k^\-^-^/'^5-^'^ 

< C4(K)A-2/'^ 

since 6 <1 and k>2. 

Define X-^- = max(0,x). Using the inequality 2{t — s)^ >t^ — and (10) 
we obtain 

P(|C - lECl >t)< 2e-2(*-|iEC-^(CI^)l)y(™.^^) + msX-' 

< c5(i^)e-^''*'(^''^^)"'"'''^'+'"'^W"' + 16VdS{K)X-^+' 
with Cd = 2-43-2'^d-'^^i/2 foj, 4.Vd{kX^^ In A)i/'^ < e{K) and any k>2. 
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